Optical manipulation of magnetism in spin-charge coupled correlated electron system 
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Photoirradiation effects in correlated electrons coupled with localized spins are studied based on 
the extended double-exchange model. In particular, we examine melting of an antiferromagnetic 
(AFM) charge order insulating state by varying the light intensity. When intense light is irradiated, 
the AFM insulating characteristics are strengthened, rather than change into the ferromagnetic 
metallic characteristic, which are expected from the conventional double exchange interaction when 
carriers are introduced by weak light irradiation or chemical doping. This provides a new principle 
for optically manipulating magnetism. 
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Optical manipulation of magnetism is one of the recent 
attractive themes in not only condensed matter physics 
but also modern nano-technologics. Since electron spins 
do not couple directly to the electric-field component of 
light, optical control of magnetism has been done in usual 
by utilizing the circular polarized light through the spin- 
orbit coupling. This strategy is applied successfully to 
the controls of magnetism in metallic magnets and some 
magnetic semiconductors. Another route for the optical 
control of magnetism is recently developed in correlated 
electron systems, in which magnetism is strongly coupled 
with electric transport. There is a possibility to change a 
magnetic state by modifying the charge state which can 
be directly controlled by light. One of the advantageous 
in this manner is a possibility of the ultrafast control 
of magnetism, since motion of electronic charge is much 
faster than pure spin motion which is usually governed 
by the exchange interaction. 

Such situation of strong correlation between mag- 
netism and electric transport is realized in a system where 
local magnetic moments arc embedded in a metallic solid. 
Through the Hund's rule coupling between the local 
moments and the conduction electrons, the ferromag- 
netic (FM) double- exchange (DE) interaction operates 
between the spatially separated local moments [HQ. Ma- 
terials in which the DE interaction gives rise to metallic 
ferromagnetism include magnetic semiconductors 0, H| , 
colossal-magnetoresistive (CMR) manganites [D, @, and 
so on. The essential aspect of the DE interaction is the 
kinetic-energy gain of conducting carriers; when local 
magnetic moments are parallel, electrons can hop be- 
tween sites without reducing the Hund's rule coupling 
energy. 

A FM metal originating from the DE interaction of- 
ten competes with other electronic phases. In particular, 
competition between a FM metal and a charge-ordered 
(CO) antiferromagnetic (AFM) insulator is widely ac- 
knowledged to be the source of the CMR effects in per- 
ovskite manganites A CO AFM insulator is charac- 
terized by an energy gap due to the Coulomb interaction 
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FIG. 1: (color online) (a): Schematic diagram of changes 
in the electronic density-of-state (DOS) and magnetic struc- 
ture on weak and strong photon excitations. Bold arrows 
represent localized spins. Broken arrows imply electronic ex- 
citations due to photoirradiation. (b): Changes in the band 
width Act 2 (bold lines) and the absolute value of the kinetic 
energy Etin (broken lines) as functions of time. Blue and red 
curves are the results for weak pumping (A p h = 0.08) and 
strong pumping (A p h = 0.35), respectively. A sign of Ekin 
is negative, and increasing (decreasing) of |-E7fcj n | implies that 
carriers gain (loss) the kinetic energy. 



and an alternating spin alignment by the AFM exchange 
interaction. Hole doping by chemical substitution of a 
CO AFM insulator induces metallic conduction with a 
FM spin alignment. Many recent experimental and theo- 
retical studies [7 - 11 have revealed that these correlated 
changes in magnetism and conduction are also realized 
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by photo-doping that introduces carriers by optical laser 
pulses, as shown in the left of Fig. [TJa). This is a naive 
extension of the DE interaction to excited states, since 
photon pumping has been widely recognized to convert 
insulators to metals in a number of organic salts 12 , l]| , 
transition- metal oxides [13, EH and other correlated elec- 
tron systems Til- 19j. 

In this Letter, we demonstrate that this conventional 
mechanism for the DE interaction does not operate in 
highly photoexcited states. As shown in Fig.QJb), when 
intense light irradiation introduces carriers to an AFM in- 
sulator, the AFM insulating characteristics are strength- 
ened, rather than change into the FM metallic character- 
istic, which are expected when it is irradiated by weak 
light. This provides a new principle for optically manip- 
ulating magnetism and conducting properties based on 
varying the light intensity. 

To examine changes in magnetism and conduction due 
to intense light irradiation, we adopt the extended DE 
model in which mobile electrons couples with local spins 
through FM Hund's rule coupling J#. The explicit 
Hamiltonian is given by 

Ho = -tJ2 + H - c -) -JhJ2 s *- S * 

■ ./ • .S' ; • \'Y."'"r ( J ) 

(ij) (ij) 



where a a is the annihilation operator for an electron 
at site i and spin cr(=t,4-), s t (= 5Z SjS ' 4s°WCi s ') and 

n i{= J2a n i° = So- c \ a c iiy) are respectively the spin and 
charge operators for electrons, and Si is the local spin 
operator. We restrict the present case to quarter-filling 
so that the average number of mobile-electrons per site is 
0.5. The first two terms correspond to the standard DE 
Hamiltonian which favors a FM metal. The competing 
phase, a CO AFM insulator, is stabilized by the Coulomb 
interaction between nearest neighbor (NN) sites, V, and 
the AFM exchange interaction between NN local spins, 
,/. In the limit of strong Hund's rule coupling, Jh 3> t, 
an electron spin and a local spin at the same site are al- 
ways parallel. This is performed by introducing the local 
SU(2) gauge transformation, and spin- less fermions, Cj. 
Then, the original Hamiltonian is reduced to the follow- 
ing spin-less fcrmion model: 
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UiCiCj + H.c. 



Vriirij 



J Si ' S i 
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where rii = cjci is the number of spin-less 
fermions, and the effective hopping amplitude 
for fermions depends on the angle of the local 
spins as %j = t{e i W i -^)/ 2 cos(^ i /2)cos(6'j/2) + 
e -i(^,-0 J )/2 s i n (g)./2) sin(6»j/2)}, as known in Ref . % 

The pump photon is introduced into the fermion 
hopping amplitude as the Peierls phase as tij —r 



tije A ( T )- dr where A(t) is the vector potential for 
the pump photon at time r, and X% IS S> position of a site 
i. The damped oscillator form is assumed for A(t) as 
A(t) = A p hee~ loT costo p hT with an amplitude A p h, a 
damping factor 70 = O.li, a photon frequency u) p h, and 
a unit polarization vector e. The time origin is taken 
to be the center of the photon pulse. The pump photon 
intensity is expressed by a dimensionless vector-potential 
amplitude A p h in terms of the light velocity, the Planck 
constant and the lattice constant. Impact by the photon 
irradiation is measured by an "effective" photon density, 
n ph = [E(t > 0) - E(t < 0)]/(Nu ph ), where E(t) is 
the electronic energy at time r, and N is the total num- 
ber of sites. The present values for the effective density, 
< riph < 0.15 for < A p h < 0.35, are realistic in the 
usual optical pump- probe experiments (2pj . 

The initial-state electronic structure, time evolution, 
and excitation spectra are calculated by the numeri- 
cal exact-diagonalization method in finite size clusters. 
The one-dimensional N(< 25) site chains and the two- 
dimensional N(= 3x5 and 5x5) site clusters with 
the open-boundary condition arc adopted. The mo- 
bile electron number is taken to be N e i e = (N + l)/2 
which corresponds to the quarter filling in the case of 
the open boundary condition. The time evolution for 
the electronic wave function, ^(r)), is obtained by solv- 
ing numerically the time-dependent Schrodinger equation 

jdMrT) 



H\i/j(t)), based on the Lanczos method [2l[. 
As for the local spin sector, we introduce the classi- 



dr 



J_pc 
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cal equations of motion given by 
d%W~- ~ §§~ = 0' where £ is the Lagrangian defined by 
C = SJ2i C0S ^i4>i ~ (?0 [22J, and (• ■ • } represents the 



average in terms of the wave function at time r. For sim- 
plicity, magnitude of spin is assumed to be S — 1/2. The 
equations of motion are solved numerically by utilizing 
the fourth order Runge-Kutta method. 

The ground-state phase diagram almost reproduces 
the previous results where the extended DE Hamiltonian 
is analyzed by the density-matrix renormalization-group 
(DMRG) method [H,!^. We chose, as an initial state of 
photo-irradiation processes, the spin-canted CO insulat- 
ing phase which is located around a boundary between 
the CO-AFM insulating and FM metallic phases. The 
photon frequency uj p h is tuned to the insulating gap en- 
ergy (see FiglU). 

Real-time spin- and charge-dynamics after photon- 
pumping were respectively monitored by the correla- 
tion function between the NN localized spins, Ks — 
N B ls ~ 2 Y,<ij) S i ■ S h and tlmt between NN mobile 
charges, Kn = N^ 1 J2{ij) (^i^j)i where Nb is the to- 
tal number of NN bonds. Time evolutions of Ksij) 
and Kn{t) for several photon amplitudes calculated in 
a TV = 15 chain are shown in Figs. a) and [UJb). The 
elapsed time r is denoted in units of t ; it is a few 
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FIG. 2: (color online): Time profiles of the NN charge 
correlation function Kn, and those of NN spin correlation 
function Ks for several pump-photon amplitudes. Results 
for the N = 15 chain cluster where V/t=lA and J/t=1.24 
are shown in (a) and (b). The pump-photon profile with 
ui p h = Lit is also shown, and photon amplitudes are A p h — 
0.04, 0.08, 0.2, 0.35. Results for the N = 25 chain (solid 
lines: V/t = 3.1, J/t = 0.625) and N = 25 square (broken 
lines: V/t = 2.1, J/t = 0.31) clusters are shown in (c) and 
(d). Pump photon energies are chosen to be 2.7t and 6.2t 
for the chain and square clusters, respectively. Shaded areas 
represent the time when the pump photons are irradiated. 



femto-seconds for conventional transition-metal oxides. 
In an ideal CO-AFM insulator, we have Kn = and 
Ks = — 1- The increase in Kn with A p h implies that the 
initial coherent CO collapses on irradiation. In contrast, 
the time profiles for Ks are rather complex. For weak 
pumping {Aph < 0.08), the AFM correlations gradually 
weaken, and exhibit oscillating behavior. This oscillation 
might be attributed to the Rabi oscillation between the 
photoexcited states where energy separations are gov- 
erned by the effective transfer, ty. These profiles in Ks 
are well understood from the conventional DE mecha- 
nism [lol . [25j , which is often applied to chemically-doped 
manganites and weakly photodopcd manganites. In con- 
trast, for strong pumping (0.2 < A p h), the AFM correla- 
tions suddenly weaken after photon pumping, and grad- 
ually recover to the initial AFM state. After rt ~ 35, the 
AFM correlation is slightly stronger than its initial value. 
In other words, for Ks when rt > 35, fewer photons col- 
lapse the AFM correlation more effectively. These char- 
acteristic features in Ks{t) and Kn{t) presented are ro- 
bust by changing the clusters, as shown in Figs.[2fc) and 
[2]Jd). Hereafter, we focus on the dynamics in the N = 15 
chain cluster to avoid redundancy and complexity. 

To investigate change in the carrier motion predicted 
by the conventional DE mechanism, Fig. [IJb) shows time 
profiles for the kinetic energy of the mobile carriers, 
E kin = (-Y^mtijCiCj + H.c.) (see broken lines). For 
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FIG. 3: (color online) Contour plots of one-particle excitation 
spectra in the momentum-energy space, (a) Results for the 
initial CO-AFM insulator at rt <C 0. White arrows denote the 
top and bottom of the valence band, (b)-(c) The spectra after 
weak photon pumping {A p h = 0.08), and (d)-(e) the spectra 
after strong photon pumping {A p h = 0.35) at rt = 5 and 35. 
White broken lines represent the top of the valence band and 
the bottom of the conduction band in the initial insulating 
state. We note that transient spectra in the region of \cj — 
1.4t| < t~ are not meaningful because of the uncertainly 
principle between time and energy. 



weak pumping {Aph = 0.08), |.E*in| increases accompa- 
nied with oscillations, as expected from the DE mech- 
anism. Completely different behavior is observed for 
strong pumping {Aph = 0.35): decreases mono- 

tonically after pumping. These results indicate that this 
highly excited state is a charge-disordered AFM insula- 
tor. Since the total energy is conserved in the simulation, 
the kinetic-energy gain for weak pumping is compensated 
by a loss in the AFM exchange-energy of the local spins. 
In contrast, for strong pumping, reductions in the ki- 
netic and Coulomb energies are compensated by gain in 
the AFM exchange-energy due to reinforcement of the 
AFM correlation. 

Time- and momentum-resolved electronic struc- 
tures provide insights into what triggers these 
unusual highly photoexcited states. Transient cor- 
related electronic bands are monitored by time- 
dependent one-particle excitation spectra, i.e. time- 
resolved (inverse) photoemission spectra (PES). 
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This is divided by the electronic part defined by 
E rn + E n + irf)~ l , and the hole part by A% u (t) = 

7' r " llm E rera l (4)r.m(4)mrfl(T)(w - E m + E t + 

ir;)" 1 , where we introduce the energy eigen value E m , the 
cigen state \m), (c k ) nm = (n\c k \m), ip n (r) = {n\ip(r)), 
and an infinitesimal constant 77 = O.lt. These formulae 
are derived based on the linear response theory in the 
photo-excited transient state 
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Results for weak 
and strong pumping are presented in Fig. [3] In the 
initial CO-AFM insulator (see Fig. [3ja)), a direct gap 
in the electron energy bands opens at the center of the 
Brillouin zone, e.g. k = n/2. After photon-pumping, 
a large change in the electronic spectra occurs due to 
many-body effects. For weak pumping at A p h = 0.08 
(left panels in Fig. [3]), changes in the electronic structure 
arc triggered by the spectral-weight transfer around 
k = 7r/2 from the valence band to the in-gap region. 
Specifically, electrons at the top of the valence band 
are excited, and a new in-gap state is generated around 
the same momentum. On the other hand, for strong 
pumping at A p h = 0.35 (right panels in Fig. [3j, the 
electronic-structure changes across the Brillouin zone. In 
particular, high spectral weights around the bottom of 
the valence band (k ~ 0) in the initial state are smeared 
out, and arc transferred to the initial gap region and the 
whole valence-band region. 

To analyze the transient valence-band width in more 
detail, we introduce the second moment of the valence 
band defined by a 2 (r) = iV" 1 J^J E fc A^Jr)^ - 
LOc) 2 duj, where uj c is the center of the band given by 
^ = I^EkAlJr^/J^^AlJr)^ and 
E top is the energy at the top of the valence band for 
rt <C 0. We define the difference from the initial band 
width as Act 2 = cr 2 (r) — cr 2 (r <C 0). Numerical re- 
sults in one-dimensional 15 site chain are presented in 
Fig. QIb). For weak pumping (A p h = 0.08), Act 2 in- 
creases and shows positive values, as expected from the 
conventional DE scenario. On the other hand, results are 
completely different for strong pumping (Aph = 0.35). 
That is, the band width becomes wide (narrow) for weak 
(strong) pumping. These are the similar behaviors with 
the changes in the kinetic energy E^m, explained previ- 
ously. 

We now present a physical picture for the highly pho- 
toexcited state (see Fig. QJa)). For weak pumping, as 
shown in Figs. [3jb) and|3jc), holes are generated near 
the top of the valence band. This situation is similar to 
the state induced by chemical doping of holes. According 
to the conventional DE mechanism, a broadening of the 
valence band accompanied with the FM alignment of lo- 
cal spins occurs in order to gain the hole kinetic-energy. 
However, for strong pumping, as shown in Figs.EJd) and 
[31(e), many photodoped holes are introduced near the 
bottom of the band. The band narrowing increases the 
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FIG. 4: (color online): Photoabsorption spectra after weak 
photon pumping (A p h = 0.08, a blue line) and strong photon 
pumping (Aph = 0.35, a red line) at time rt — 35. Absorption 
spectra for the initial CO-AFM insulator are also shown by 
a black line. Bold arrow indicates the pump photon energy 
ujph = Lit which is a peak energy of the optical absorption 
spectra in the initial state. We note that transient spectra in 
the region of |oj < 0.03i are not meaningful because of the 
uncertainly principle between time and energy. 



hole kinetic-energy gain. This is considered to be the ori- 
gin of the enhanced AF insulator characteristics induced 
by strong photon pumping. 

We propose that time-resolved pump-probe experi- 
ments can confirm the present unusual photoexcited 
state. The pump-probe optical absorption spectra is 
calculated by the linear response theory for the probe 
photon defined by a(ui; r) = — 7r^ 1 Im(J r (a; — H + (TL) + 
where we introduce the current operator given 
by J = Ujcfej + H.c. [Hj]. Figure H shows op- 

tical pump-probe spectra (the transient optical absorp- 
tion spectra a(w;r)) at time rt = 35. A gap-like fea- 
ture is observed in a(u>; r) for an initial CO insulator. 
A low-energy peak grows in the initial optical gap for 
Apt = 0.08, as expected for the metallic state. On the 
other hand, the reduction in the whole spectral inten- 
sity in Aph = 0.35 implies enhancement of the localized 
character of mobile carriers. The low spectral intensity 
in the initial gapped region for strong pumping indicates 
collapse of the coherent CO state. 

The present study reveals a new aspect of the DE sys- 
tem. The reinforcement of the AFM correlation and the 
band narrowing observed in the present study were not 
expected for chemical doping and thermal carrier cre- 
ation. These phenomena are only realized in photoirra- 
diated spin-charge coupled systems for which high elec- 
tronic excitation is possible at high photon densities, and 
the band width is directly controlled by the spin config- 
urations. Thus, the strongly pumped state is identified 
as a photo-induced hidden state, which will be directly 
detected by recently developed time-resolved resonant 
x-ray diffraction and angular resolved PES. The clear 
difference between weak and strong pumping provides 
a possible means for optically manipulating magnetism 



5 



and conduction based on photon intensity. After several- 
hundred femto seconds, which is beyond the limits of 
the present simulation, spin relaxation originating from 
the spin-orbit interaction begins to occur and it tends to 
stabilize the FM metal. Energy relaxation will eventu- 
ally restore the system to the initial CO-AFM insulating 
state. In other words, relaxation effects for both angular 
momentum and energy relax the system to trivial fixed 
points. The unique photo-induced phenomena occur in 
early stages before relaxation starts to operate. 
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